Field-induced criticality in a gapped quantum magnet with bond disorder. 
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Neutron diffraction and calorimetric measurements are used to study the field-induced quantum 
phase transition in piperazinium-Cu2(Cli_^Br I )6 (x = 0, x = 3.5% and x = 7.5%), a prototypical 
quantum antiferromagnet with random bonds. The critical indexes <j> and /3 are determined. The 
findings contradict some original predictions for Bose Glass based on the assumption z = d, but are 
consistent with recent theoretical results implying z < d. Inelastic neutron experiments reveal that 
disorder has a profound effect on the lowest-energy magnetic gap excitation in the system. 

PACS numbers: 75.10.Jm, 64.70.Tg, 75.40.Cx, 72.15.Rn 



The persistent interest in quantum magnets is driven 
by their use as prototypes for the study of quantum criti- 
cal phenomena. In this context, much attention has been 
given to field-induced phase transitions in gapped spin 
systems with a singlet ground state, the so-called quan- 
tum spin liquids. Such transitions are described in terms 
of a Bose Einstein condensation (BEC) of magnonsPThe 
specific advantage of magnetic realizations of BEC lies in 
the relatively easy experimental access to- and control of- 
the key relevant quantities, notably the effective chemical 
potential, boson density and BEC order parameter. In 
magnetic BEC, the latter correspond to magnetic field, 
magnetization and antiferromagnetic order, respectively. 

A particular area of recent interest is BEC in the pres- 
ence of a spatially random potential. Such disorder leads 
to a localization of bosons and the appearance of a novel 
quantum phase, the so-called Bose GlassPS* Most im- 
portantly, it changes the universality class of the quan- 
tum critical point. While in the disorder-free case it is 
Mean Field (MF) like with dynamical critical exponent 
z = 2P very unusual scaling with z equal to dimension- 
ality (z = d) was predicted for disordered systemsPThis 
premise was recently challenged, 5 6 suggesting that in 
fact z < Quantum magnetic materials offer a unique 
opportunity of addressing this important matter of de- 
bate experimentally. A random potential for magnons 
can be created by randomizing the strength of magnetic 
interactions. This is achieved through a spatially ran- 
dom chemical modification. Indications of the magnetic 
Bose GlassSHDJ. anc j a modification of the field-induced or- 
dering transitiorP^HS] have been reported. However, ex- 
periments only added to the confusion regarding critical 
indexes and scaling. Recent studies of TI0.64K0.36CUCI3 
seemed to validate the z = d hypothesis.^ Other ex- 
periments, on the organic compound DTN,^ yielded 
quite different critical exponents, and were supported 
by Density Matrix Renormalization Group (DMRG) cal- 
culations. In the present work we seek to resolve the 



issue through neutron and calorimetric studies of the 
chemically disordered quantum magnets piperazinium- 
Cu2(Cli_ 3; Br a ;)e (PHCX). We determine the critical ex- 
ponents <f> an d P, the latter previously unmeasured for 
this type of transition. We find experimental values that 
are clearly inconsistent with predictions based on z = d. 

One rationale for selecting PHCX as a test material 
is that the disorder-free Cl-rich (x = 0) compound, 
(PHCC) is an exception ally well-characterized Heisen- 
berg quantum magnet E2HI2I It features a complex quasi 
two-dimensional spin network of S = 1/2 Cu 2+ ions 
bridged by Cu-Cl-Cl-Cu superexchange pathways, see 
Fig.l in ReffT5l The exact Hamiltonian for this sys- 
tem remains unknown with 6 to 8 possible superexchange 
pathways within the CU2CI6 layers. However, for the 
physics discussed here relevant ramification is the spin 
liquid ground state and the gapped (A = 1 meV) magnon 
spectrum with bandwidth of about 1.7 meV. In magnetic 
fields exceeding noH c — 7.5 T, PHCC undergoes a quan- 
tum phase transition to a state with spontaneous anti- 
ferromagnetic (AF) order perpendicular to the direction 
of applied field. The critical properties of this transition 
have been investigated using a variety of techniques.^ 
The two critical exponents studied in this work are the 
order parameter exponent /? and the so-called crossover 
exponent <f/. The former defines the field dependence of 
the AF ordered moment at T 0: |(Sj_)| ~ (H - H C )P, 
where H c = H C (T = 0). The crossover index <fi defines 
the phase boundary of the ordered state on the H — T 
phase diagram: T C {H) ~ (H - H c )^ or T^H) 1 ^ ~ 
(H-H c ). 

Br substitution in PHCX is expected to locally af- 
fect the bond angles in the halogen-mediated superex- 
change routes, and thereby to realize the random bond 
model. This mechanism was previously exp loite d in 
the study of random-bond systems IPA-CuClaJ 211 ^ Sul- 
Cu 2 Cl 4 (Ref. ED and DTNP High quality PHCX sin- 
gle crystals were grown from solution as described in 
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Ref. UHJ Single crystal X-ray diffraction confirmed that 
for Br content x up to at least 10%, the crystal symmetry 
remains intact, space group (PI). The lattice parameters 
are a = 7.9691(5)1,6 = 7.0348(5)A,c = 6.0836(4)A,a = 
111.083(3)°, = 99.947(3)°, 7 = 81.287(4)° for the pure 
material with x = 0p^ and the cell dimensions increase 
linearly with x. Here x stands for nominal Br concen- 
tration in the starting solution. Studies on similar com- 
pounds have shown a very good correspon dence between 
nominal and final dopant concentrations.^^ Most of 
the work reported here was done on perfectly faceted 
single crystals with x = 3.5 % and x — 7.5 %, where 
a = 7.9799(2) A,b = 7.0406(2)A, c = 6.09090(10) Aand 
a = 7.9905(3) A, b = 7.0454(2)A, c = 6.1000(2) A, re- 
spectively, with the triclinic angles indistinguishable from 
the pure compound's. 

The crossover exponent was derived from calorimctric 
measurements. For x = 0, x — 3.5 % and x — 7.5 %, the 
data were collected using a 14 T Quantum Design PPMS 
with a 3 He- 4 He dilution insert, both as a function of tem- 
perature and magnetic field applied along the o* axis. 
Typical scans are shown in Fig{T^,,b (symbols). Sharp 
lambda-anomalies are apparent in all cases and mark the 
ordering transition. The apparent broadening at low T 
for x = 7.5 % sample is of instrumental origin. It is 
caused by differences in the slope of the phase boundary 
in conjunction with finite temperature resolution. In the 
constant-field scans, the actual transition temperatures 
T c were determined by fitting the data at each field, in 
the vicinity of the peaks, to an empirical form of power 
law with broadening!-^ A similar approach was used to 
obtain the temperature dependence of the critical field 
H c from constant-temperature scans. In this case the fit 
function was a composite of two half-Lorentzians of non- 
equal widths and common maxima. The resulting phase 
boundaries are plotted in Fig. [2] (symbols). Clearly, the 
introduction of disorder shifts the transition to higher 
fields and lower temperatures. 

Any more quantitative analysis has to take into ac- 
count the experimental errors. For consistency, we only 
used the more comprehensive data set obtained in con- 
stant temperature scans (empty symbols in Fig. [2]) . The 
error bar on H c was estimated as the standard devia- 
tion of this parameter in the least squares fits described 
above and is significantly smaller than the symbol size 
in the plot. The dominant error on T c is of instrumental 
origin, and due to the use of the temperature relaxation 
measurement method. It is defined by the magnitude of 
applied heat pulse (3% of the starting temperature in our 
case). 

The exponent <p is obtained by weighted least squares 
fit of the power law function to the experimental field- 
dependence of T c accounting for experimental errors on 
both. This deceptively simple procedure is actually quite 
delicate. The fit is performed in a progressively shrinking 
data range T < T max — > 0. In Figs. [2jD-[2ji, we show the 
fitted values of <j> for the three materials plotted against 
Pmax- The range used for the final fits is to be selected 
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FIG. 1. (color online) Symbols: Heat capacity of PHCX with 
varying Br concentration x measured (a) as a function of mag- 
netic field at T=0.5K and (b) as a function of temperature at 
H = 14 T. The solid lines are fits to the data as described in 
the text. 
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TABLE I. Experimental estimates for the critical field H c in 
H||a* orientation and the critical exponents cj> (from calori- 
metric data) and j3 (from neutron diffraction) of the field- 
induced quantum phase transition in PHCX in comparison 
with recent measurements on similar systems and theoretical 
results, whose relevance is discussed in the text. 



based on two conflicting requirements. It has to be as 
narrow as possible, in order to access the true critical 
region. This desire is counteracted by the divergence of 
confidence intervals at T max — > 0. As a rule of thumb, 
one can choose the largest T max that gives a value of <J) 
within the confidence intervals obtained using any nar- 
rower ranges. Further reducing the fitting range does not 
produce a statistically significant change in the result. 

Returning to our case of PHCC, based on the plots 
in Fig. [2] and the above-mentioned rule of thumb, for the 
x = 0, x = 3.5% and x = 7.5% materials we have selected 
fitting ranges T max = 1.75 K, 0.9 K and 1.15 K, respec- 
tively. The resulting values of <fr and H c are summarized 
in Table |Ij Power laws derived from these parameter val- 
ues are shown as solid lines up to T = T max in Fig. [2j 

The behavior of the magnetic BEC order parameter 
was studied in neutron diffraction experiments. For the 
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FIG. 2. (color online) (a) Phase boundaries of ordered state 
for Br-free, 3.5% and 7.5% Br substituted PHCX. Open and 
filled circles are from C(H) and C(T) measurements, respec- 
tively. Solid lines are a power laws H — H c ~ T C {H) X '^ with 
critical fields H c and crossover exponents <f> derived from the 
data as described in the text, (b-d) Crossover exponent cf> 
obtained from least squares fit as a function of fitting window 
given by T < T max . Filled symbols show the value quoted in 
Table [TJsatisfiying our best fit criteria described in the text. 



x = 3.5% and x = 7.5% samples these were carried out at 
the CRG-CEA D23 diffractometer at ILL and the RITA- 
II 3-axis spectrometer at PSI, respectively. In the two ex- 
periments, we used Cu monochromator to produce neu- 
trons of incident energies 50 meV at ILL and pyrolitic 
graphite for 4.6 meV neutrons with Be filter in front of 
the analyzer at PSI. Sample environment was in all cases 
a split coil cryomagnet with a 3 He- 4 He dilution cryostat. 
The magnetic field was applied along a* in the first and 
along the b crystallographic axis in the second experi- 
ment. As previously reported for PHCC, 17 long range 
order in the BEC phase is marked by the appearance of 
new magnetic Bragg reflections half-integer reciprocal- 
lattice points (h + 1/2,0,/ + 1/2). The inset in Fig^ 
shows typical scans across the (0.5, 0, —0.5) position mea- 
sured in the x = 7.5% sample at T = 0.5 K in the BEC 
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FIG. 3. (color online) Upper panels: Neutron diffraction in- 
tensity of the (|0|) peak measured in PHCX, a) as a function 
of magnetic field at T=160mK, b) as a function of tempera- 
ture. Open and solid symbols correspond to zero field cooling 
and field cooling, respectively. Squares (circles) correspond 
to 3.5% (x — 7.5%) Br content. Lower left corner: the trajec- 
tory in (H,T) space for the measurements on the x = 7.5% 
material. Inset: rocking curve measured at T = 0.5 K (filled) 
and T = 2 K (empty circles) in that sample. Lower panels: 
order parameter exponents as obtained by shrinking the fit- 
ting interval analysis of the 3.5% (panel c) and 7.5% Br (panel 
d) samples data. Filled symbols show the best fit values of j3. 



phase, and the background collected at T — 2 K in the 
paramagnetic state. To within experimental error, all 
magnetic reflections are resolution-limited. 

The critical index /3 was determined from measure- 
ments of the (0.5, 0, —0.5) Bragg peak intensity in zero 
field cooled (ZFC) samples. These data were taken upon 
increasing the applied field at T = 160 mK and are plot- 
ted in Figj3ji. /? was then determined in power law fits, 
assuming that peak intensity scales as as the square of the 
order parameter. As for (f>, the fits for j3 were performed 
for a series of shrinking fitting intervals, of magnetic field 
SH = H — H c in this case. The fit results and standard 
deviation confidence intervals are plotted against SH in 
Figs. |3^,d. Following the reasoning described earlier, in 
our final analysis we have selected ^qSH = 1.86 T and 
2.15 T, for x = 3.5% and x = 7.5%, respectively. The 
resulting exponents are summarized in Table [I] In order 
to properly characterize the quantum critical point, one 
has to know j3 in the limit T — >• 0. Previous studies of 
PHCC suggest that f3 remains practically constant for 
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E, meV 

FIG. 4. (color online) Inelastic neutron scattering intensity at 
(0.5,0.5,-0.5) in energy transfer range from 0.5 to 2.5 meV for 
3.5% (green squares) and x — 7.5% Br (red circles) content 
samples. Corresponding solid lines are the fitted Lorentzian 
energy profiles convoluted with the 4-dimensional instrumen- 
tal resolution function. For 0.0% Br sample resolution con- 
voluted simulation is shown for comparison using parameters 
from ReffT5l 

T < 0.3 KP^ We have additionally verified this behavior 
for our 3.5% PHCX sample. Thus, the 160 mK measure- 
ments are good estimates of the zero-temperature values. 
Comparing our results with those reported in Ref. 1171 we 
conclude that disorder does not significantly change the 
criticality of the order parameter. 

This is not to say that the effect of disorder is negligi- 
ble. Inelastic neutron scattering experiments performed 
at zero field on the TASP 3-axis spectrometer^! at PSI 
using 3.5meV incident-energy neutrons and Be filter re- 
veal that while the magnons are sharp in PHCC, 113 they 
acquire a considerable intrinsic width in disordered sam- 
ples. Typical constant-q scans at the AF zone-center 
(0.5,0.5,-0.5) in the x=3.5% and x=7.5% samples are 
shown in Fig. [3] by the open symbols. The solid lines 
are fits of Lorentzian energy profiles convoluted with the 
4-dimensional instrumental resolution functional For 
x=0% the line is a simulation assuming that peak shape 
is determined solely by resolution. The observed increase 
of gap energy with Br substitution is fully consistent 



with the observed increase of critical field and may be 
attributed to a change in the disorder-averaged exchange 
constants. However, the large intrinsic energy width 
T = 0.1 meV and 0.25 meV for x=3.5% and x=7.5% sam- 
ples, respectively, is entirely a disorder-induced effect!^. 

Our present diffraction experiments show that disor- 
der also affects the behavior of the order parameter deep 
within the high-field phase. In the 7.5% material the 
intensity of the magnetic reflection exhibits a history 
dependence. This is illustrated in Fig. |3^,,b. Here we 
show data collected along a trajectory illustrated in the 
inset. No such behavior was reported for the disorder- 
free material, but a similar, and much stronger history 
effect was previously studied in a related system IPA- 
(Cui^Br^ClaPS Related phase transition scenario that 
could lead to hysteretic behavior was discussed in Ref. [521 

Theoretical and numerical predictions for the two mea- 
sured critical exponents are shown in Table [I] for a direct 
comparison with our experiments. In the absence of dis- 
order, one expects MF values^ (3 = 1/2 and 4> — 2/3, in 
reasonable agreement with the data for PHCC. Assum- 
ing z = d the disordered case gives much larger value 
for the crossover exponent (f> — 2p^ This prediction is 
clearly inconsi stent with our experiments. More recent 
numerical worlP^^ suggested smaller values of the ex- 
ponent. Those from Refs. [7] and [5] are less relevant to 
3-dimcnsional ordering in PHCX, as they are obtained 
for a 2-dimensional model. On the other hand, the QMC 
calculations of Ref. 13 should be directly comparable to 
our data. The corresponding estimates of 4> are m a 
rather good agreement with our experiments, as well as 
with those on the Br-substituted DTN compound!^ A 
relevant theoretical study of /3 is presently lacking. The 
claimed observation^ of values consistent with z = d can 
be due to overinterpretation of limited dataset. Evidence 
is mounting for the breakdown of the z = d hypothesis 
in the Bose Glass to BEC quantum phase transition. 
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